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43 ■ ABSTRACT 

a: 

Inelastic photoproduction of J/ip particles at high energies is one of the processes to 
determine the gluon distribution in the nucleon. We have calculated the QCD radiative 
corrections to the color-singlet model of this reaction. They are large at moderate photon 
^ energies, but decrease with increasing energies. The cross section and the J/ip energy spec- 

trum are compared with the available fixed-target photoproduction data and predictions are 
given for the HERA energy range. 
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The measurement of the gluon distribution in the nucleon is one of the important goals of 
lepton-nucleon scattering experiments. The classical methods exploit the evolution of the 
nucleon structure functions with the momentum transfer, and the size of the longitudinal 
structure function. With rising energies, however, jet physics and the production of heavy 
quark states become important complementary tools. Besides open charm and bottom pro- 
duction, the formation of J/ip bound states |J in inelastic photoproduction experiments 

7 + N^J/ip + X (1) 

provides an experimentally attractive method since J/ip particles are easy to tag in the 
leptonic decay modes. 

Many channels contribute to the generation of J/ip particles in photoproduction experi- 
ments H, similarly to hadroproduction experiments. However, no satisfactory quantitative 
picture has emerged yet and the production of a large surplus of ip' particles in pp collisions 
awaits the proper understanding. Theoretical interest so far has focussed on two mecha- 
nisms for J/ip photo- and electroproduction, elastic/diffractive [[3],|J and inelastic production 
through photon-gluon-fusion Jl|,0. While by the first mechanism one expects to shed light 
on the physical nature of the pomeron, inelastic J/ip production provides information on the 
distribution of gluons in the nucleon [H . The two mechanisms can be separated by measuring 
the J/ip energy spectrum, described by the scaling variable 

z = ^ (2) 

with p, fc^, )7 being the momenta of the nucleon and J/ip, 7 particles, respectively. In the nu- 
cleon rest frame, z is the ratio of the J/ip to the 7 energy, z = E^/E 1 . For elastic/diffractive 
events z is close to one; the inelastic events are experimentally restricted in general to the 
range z < 0.9 ||. The production of J/ip particles at large transverse momenta is dominated 
by gluon fragmentation mechanisms JF]. 

Inelastic J/ip photoproduction through photon-gluon fusion is described in the color- 
singlet modelQ through the subprocess 

7 + g -+ J/i) + g (3) 

shown in Fig.l. Color conservation and the Landau- Yang theorem require the emission of a 
gluon in the final state. The cross section is generally calculated in the static approximation 
in which the motion of the charm quarks in the bound state is neglected. In this approxima- 
tion the production amplitude factorizes into the short distance amplitude 7 + g — > cc + g, 
with cc in the color-singlet state and zero relative velocity of the quarks, and the cc wave 
function y?(0) of the J/ip bound state at the origin which is related to the leptonic width. 
The cross section of the subprocess (^) may be written as 

dqW 128tt 2 aa y c 2 |y(0)| 2 s 2 s 2 + t 2 t 2 + u 2 u 2 
dh ~ 3 s 2 J ^ M m s\t\u\ {) 



1 Dual approaches to J/ip photoproduction based on the subprocess 7 + g — * c + c are, in contrast to their 
e~ analogon, not sharply defined and they will therefore not be discussed in the present context. 



where t denotes the momentum transfer squared from the photon to the J/ip particle, and s,u 
are the total energy and the momentum transfer from the photon to the gluon \t\ — t — M'j,, 
etc]. The cross section (ffl) is infrared finite. The photoproduction cross section on the 
nucleon is obtained by integrating the subprocess over the gluon flux, 



<r = /d£<7(0* (5) 



When confronted with photoproduction data of fixed-target experiments |||J , the theoretical 
predictions underestimate the measured cross section in general by more than a factor two, 
depending in detail on the J/ip energy and the choice of the parameters [for a recent review 
see Ref. @]. The discrepancy with cross sections extrapolated from electroproduction data 
ToUTTj is even larger. 



The lowest-order approach to the color-singlet model demands several theoretical refine- 
ments: (i) Higher-order perturbative QCD corrections; (ii) Relativistic corrections due to the 
motion of the charm quarks in the J/ip bound state; and last but not least, (iii) Higher-twist 
effects which are not strongly suppressed due to the fairly low charm-quark mass. While the 
relativistic corrections have recently been demonstrated to be under control in the inelastic 
region [jTJ], the problem of higher- twist contributions has not been approached so far. The 
analysis of the higher-order perturbative QCD corrections will be presented in this note. Ex- 
pected a priori and verified subsequently, these corrections by far dominate the relativistic 
corrections in the inelastic region, being of the order of several a s (Mj/ , ) ~ 0.3. In the first 
step of a systematic expansion, they can therefore be determined in the static approach [|13j . 

Generic diagrams which build up the cross section in next-to-leading order are depicted 
in Fig.l. Besides the usual self-energy diagrams and vertex corrections for photon and gluons 
(b), we encounter box diagrams (c), the splitting of the final-state gluon into gluon and light 
quark- ant iquark pairs, as well as diagrams renormalizing the initial-state parton densities 
(e). The evaluation of these amplitudes has been performed in the Feynman gauge. We 
have adopted the dimensional regularization scheme to calculate the singular parts of the 
amplitudes. [We have refrained from using the seemingly simpler scheme of dimensional 
reduction which however gives rise to complications in the massive quark case as pointed 
out in Ref. [0.] The masses of light quarks in Fig.l(d,ei) have been neglected while the 
mass parameter of the charm quark has been defined on-shell, with m c = ^Mj/^p for consis- 
tency. Since the proper renormalization scale lir is expected in the range of about m c to 
Mj/^p = 2m c , we have used the MS scheme ||15P with four active flavors for the renormaliza- 



tion of the QCD coupling a s nf {(J>r)- We have carried out the renormalization program also 
in the extended MS scheme |T(| in which the massive particles are decoupled smoothly for 



momenta smaller than the quark mass. The results in the two schemes differ marginally, by 
typically less than 6 %. The exchange of Coulombic gluon quanta in the diagram (lc) leads to 
a Coulomb singularity ~ ir 2 /2j3r which can be isolated by introducing a small relative quark 
velocity (3r. Following the standard path fL7| |, we interpret this effect as Sommerfeld rescat- 



tering correction which can effectively be mapped into the cc wave function. As expected, 
the infrared singularities cancel when the emission of soft and collinear final-state gluons and 
light quarks, characterized by a cut-off A [|14| , |18| , is added to the virtual corrections. The 



collinear initial-state singularities can be absorbed, as usual, into the renormalization of the 



parton densities [19] defined in the MS factorization scheme. 



After carrying out this straightforward but tedious program, the results for the cross 
sections of the subprocesses can be presented in the form of scaling functions, 
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i = 9^-iQ denoting the parton targets. For the sake of simplicity, we have identified the 
renormalization scale with the factorization scale p? R = fi 2 F = Q 2 . The scaling functions 
depend on the energy variable r\ = s/Am 2 — 1. c^J is the lowest-order contribution which 
scales ~ i]^ 1 ~ Am 2 / s asymptotically c^ can be decomposed into a "virtual + soft" (V+S) 



piece, and a "hard" (H) gluon-radiation piece; the In- 7 A singularities of the (V+S) cross 
section are mapped into (H), cancelling the equivalent logarithms in this contribution so 
that the limit A — > can safely be carried out [|H|. The nomenclatura "hard" and "virtual 
+ soft" is therefore a matter of definition, and negative values of c^ may occur in some 
regions of the parameter space. Up to this order, the wave-function at the origin is related 
to the leptonic J/ip width by 
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with only transverse gluon corrections taken into account explicitly 



The scaling functions c 7 i(rj) are shown in Figs.2a/b for the parton cross sections integrated 
over z < z\ where we have chosen z\ = 0.9 as discussed before. [Note that the definition of 
z is the same at the nucleon and parton level since the momentum fraction £ of the partons 
cancels in the ratio (§).] The following comments can be inferred from the figures, (i) The 
form of the scaling functions resembles the scaling functions in open-charm photoproduction 
||18|| . However, there is an important difference. The "virtual + soft" contribution for J/ip 
production is significantly more negative than for open-charm production. The destructive 
interference with the lowest-order amplitude is not unplausible though, as the momentum 
transfer of virtual gluons has a larger chance [in a quasi-classical approach] to scatter quarks 
out of the small phase-space element centered at p c + p^ = pj/^ than to scatter them from 
outside into this small element, (ii) While c^J and c^ ,v+s ) scale asymptotically ~ 1/s, the 
hard coefficients c^' H ) and c^ [as well as c^] approach plateaus for high energies, built-up 
by the flavor excitation mechanism, (iii) The cross sections on the quark targets are more 
than one order of magnitude smaller than those on the gluon target, (iv) A more detailed 
presentation of the spectra would reveal that the perturbative analysis is not under proper 
control in the limit z — > 1, as anticipated for this singular boundary region [21]. Outside the 
diffractive region, i.e. in the truly inelastic domain, the perturbation theory is well-behaved 
however. 

The cross sections for J/ip photoproduction on nucleons are presented in Figs. 3a and 
3b. In the first of the two figures we compare the leading-order and next-to-leading order 
calculations with the J/ip energy spectra of the two fixed-target photoproduction experiments 
at photon energies near EL = 100 GeV, corresponding to invariant energies of about y/s 



ip 



14 GeV. The GRV parametrizations of the parton densities |22[ have been used. Since the 
average momentum fraction of the partons < £ >~ 0.1 is moderate, the curves are not 
sensitive to the parametrization in the small-x region. The results are shown for two values 



of a^(Mj,^) = 0.25 and 0.30 which correspond to the average fit value in Ref.|^3[ and 
the la upper boundary of the error band, respectively. Since the cross section depends 
strongly on the QCD coupling, we adopt this measured value, thus allowing for a slight 
inconsistency to the extent that the GRV fits are based on a marginally lower value of a s . 
The ^-factor K = onlo/^lo turns out to be ~ 2.45 with one part ~ 1.73 due to the QCD 
radiative corrections of the leptonic J/ip width and a second part ~ 1.44 due to the dynamical 
QCD corrections per se. The if -factor is nearly independent of z. The dependence on the 
renormalization/factorization scale Q is reduced considerably in next-to- leading order pTf . 
While the ratio of the cross sections in leading order for Q = m c : (V2 m c ) : Mj/^ is given by 
1.7 : 1.3 : 1, it is much closer to unity, 0.9 : 1.1 : 1, in the next-to-leading order calculation, 
Fig. 4. The overall behavior of the Q dependence is reminiscent of heavy quark production in 
hadron collisions P4| . The cross section runs through a maximum |25[ near Q « 1.4 m c with 



broad width, the origin of the stable behavior in Q. In the BLM scheme |26j Q moves from 
values below m c at low energies up to ~ 1.5 m c at the HERA energy of y/s 1P ~ 100 GeV. 
In particular the value at high energies is significantly larger than the corresponding BLM 
value for J/ip decays. The typical kinematical energy scale is not set any more by the small 
gluon energy in the J/ip decay but rather by the typical initial-state parton energies. 

In a systematic expansion we may finally add the relativistic corrections as estimated 
in Ref.|[12||. Two conclusions can be drawn from the final results presented in Fig. 3a. (i) 



The J/ip energy dependence da/dzi^ + N — > J'/ip + X) is adequately accounted for by 
the color-singlet model so that the shape of the gluon distribution in the nucleon can be 
extracted from J/ip photoproduction data with confidence, (ii) The absolute normalization 
of the cross section is somewhat less certain; this is apparent from the comparison with 
the photoproduction data. [The situation is worse for electroproduction data [0]. However, 
allowing for higher-twist uncertainties of order (A/m c ) k <^ 20% for k > 1, we conclude that 
the normalization too appears to be under semi- quantitative control - at the least. 

In Fig. 3b we present the prediction of the cross section for the HERA energy range. In 
this high energy range the K-iactor is smaller than at low energies, a consequence of the 
negative dip in the c^ scaling function of Fig. 2. For a^(Mj, ,) = 0.25 we find, for z < 0.9, 
a value of about a{^+p — > J/ip+X) w 21 nb at an invariant 7p energy of y/s 1P ~ 100 GeV; 
this value rises to about 32 nb if we choose the larger value 0.30 for the QCD coupling. The 
production of T bottonium bound states is suppressed by a factor of about 300 at HERA, 
a consequence of the smaller bottom electric charge and the phase space reduction by the 
large h mass. 

In summa. We have shown in this next-to-leading order perturbative QCD analysis that 
the energy shape of the cross section for J/if; photoproduction is adequately described by 
the color-singlet model. A semi- quantitative understanding, at the least, has been achieved 
for the absolute normalization of the cross section. Higher-twist effects must be included to 
improve the quality of the theoretical analysis further. The predictions for the HERA energy 
range provide a crucial test for the underlying picture as developed so far in the perturbative 
QCD sector. 
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Figures 

Fig. 1. Generic diagrams for inelastic J/ip photoproduction: (a) leading order contribution; 
(b) vertex corrections; (c) box diagrams; (d) splitting of the final state gluon into 
gluon or light quark-antiquark pairs; (e) diagrams renormalizing the initial-state 
parton densities. 

Fig. 2. (a) Coefficients of the QCD corrected total inelastic [z < 0.9] cross section 7 + g — > 
J/ip + X in the physically relevant range of the scaling variable rj = s lp /4m 2 — 1; 
and (b) for 7 + q/q — > J/ijj + X. 

Fig. 3. (a) Energy spectrum da/dz, at the initial photon energy E^ = 100 GeV, for two 
different values of a s (Mj,^) compared with the photoproduction data [§]|§. (b) 
Total cross section for inelastic J/if; photoproduction 7 + P — > J/ip + X as a 
function of the photon-proton center of mass energy in the HERA energy range. 
[The dashed line presents the T photoproduction cross section, amplified by a factor 
100.] 

Fig. 4. Dependence of the total cross section 7 + g — > J/ip + X on the renormaliza- 
tion/factorization scale Q at a photon energy of E 1 = 100 GeV. 
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